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Abstract 

Multi-agent systems cooperation to achieve global goals is usually limited by sensing, actuation, and 
communication issues. At the local level, continuous measurement and actuation is only approximated 
by the use of digital mechanisms that measure and process information in order to compute and update 
new control input values at discrete time instants. Interaction with other agents or subsystems takes 
place, in general, through a digital communication channel with limited bandwidth where transmission 
of continuous-time signals is not possible. Additionally, communication channels may be subject to 
other imperfections such as time-varying delays. This paper considers the problem of consensus (or 
synchronization of state trajectories) of multi-agent systems that are described by general linear dynamics 
and are connected using undirected graphs. An event-triggered consensus protocol is proposed, where 
each agent implements discretized and decoupled models of the states of its neighbors. This approach 
not only avoids the need for continuous communication between agents but also provides a decentralized 
method for transmission of information in the presence of time-varying communication delays where 
each agent decides its own broadcasting time instants based only on local information. This method 
gives more flexibility for scheduling information broadcasting compared to periodic and sampled-data 
implementations. The use of discretized models by each agent allows for a periodic event-triggered 
strategy where continuous actuation and continuous measurement of the states are not necessary. 
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I. Introduction 


Cooperative control of multi-agent systems is an active research area with broad and rele¬ 
vant applications in commercial, academic and military areas [HI. The design of decentralized 
and scalable control algorithms provides the necessary coordination for a group of agents to 
outperform agents operating independently. In multiple scenarios where communication among 
agents is limited, decentralized computation of the time instants that each agent needs to transmit 
relevant information is also necessary. Continuous actuation and continuous measurement of local 
states may also be restricted by particular hardware limitations. 

Consensus problems with limited communication have been studied using the sampled-data 
(periodic) approach [|3, [[3J, 0, and [|5]]. An important drawback of periodic transmission is 
that it requires synchronization between the agents, that is, all agents need to transmit their 
information at the same time instants and, in some cases, it requires a conservative sampling 
period for worst case situations. 

In the present paper, in lieu of periodic approaches, we use an asynchronous communication 
scheme based on event-triggered control strategies and we consider agents that are described by 
general linear dynamics which are subject to limited actuation update rates and also to limited 
local sensor measurement update rates. In addition, we consider the case where communication 
among agents is subject to communication delays. In contrast to periodic (or time-triggered) 
implementations, in the context of event-triggered control, information or measurements are not 
transmitted periodically in time but they are triggered by the occurrence of certain events. In 
event-triggered broadcasting d6j, Gl, [O, [O, ifTOl . ifTTIl . Ifl2l . and |[T3ll . a subsystem sends its local 
state to the network only when it is necessary, that is, only when a measure of the local subsystem 
state error is above a specified threshold. Event-triggered control strategies have been used for 
stabilization of multiple coupled subsystems as in |[T4ll . [fl5l . and [fT6ll . Consensus problems have 
also been studied using these techniques lH7l . Ifl8l . Ifl9l , lf20l . lITTl . |[22l . If23l . Event-triggered 
control provides a more robust and efficient use of network bandwidth. Its implementation in 
multi-agent systems also provides a highly decentralized way to schedule transmission instants 
which does not require synchronization compared to periodic sampled-data approaches. 

One important restriction related to event-triggered control techniques is that continuous 
measurement of state variables and continuous computation of state errors and time-varying 
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thresholds is required. One solution explored by different researchers is self-triggered control 
[21, ll25l . Il26l . Il27l . The main difference with respect to event-triggered control is that a measure 
of the state is not being compared constantly against a predefined threshold. Instead, the current 
state measurement is used to determine its next deadline, i.e. the next time that the sensor is 
required to send a measurement to the controller. A recent extension to event-triggered control 
is the so called periodic event-triggered control ll28ll . If29l . where measurements of states and 
computations of errors and thresholds occur, not continuously in time, but only at periodic time 
instants. The decision to transmit the current measurement is taken at the sampling instants based 
on current measurements and computed thresholds. 

Consensus problems where all agents are described by general linear models have been 
considered by different authors EOt EH . Ii32ll . E3l . E4il . E51 . E6l . and E71 . In these papers 
it is assumed that continuous communication between agents is possible. The work in ll38il con¬ 
siders the consensus problem of agents with linear dynamics under communication constraints. 
Specifically, the authors consider the existence of continuous communication among agents for 
finite intervals of time and the total absence of communication among agents for other time 
intervals, and the minimum rate of continuous communication to no communication is given. 

Event-triggered consensus of agents with linear dynamics and limited communication was 
recently explored in E9l and li40l . In our previous work PTll and ll42il we proposed a novel 
approach in which each agent implements models of the decoupled dynamics of each one of its 
neighbors and uses the model states to compute the local control input. This approach offered 
better performance than Zero-Order-Hold (ZOH) approach used in E9l and BOl where the 
updates from neighbors are kept constant by the local agent. A similar model-based framework 
was proposed in PEI where only constant thresholds were used. One of the main limitations of the 
ZOH approach E9ll . flOl is that it is not capable to keep up with unstable trajectories and updates 
need to be generated more frequently. In consensus with general linear dynamics, unstable 
systems are one of the most interesting cases to analyze. The model-based approach in iHTl 
lf42l provides better estimates of neighbors and reduces generation of events as agents converge to 
similar unstable trajectories. Communication delays using the event-triggered approach in these 
papers were addressed in P4l . The present paper offers complete proofs with respect to PHI and, 
more importantly, it extends the consensus protocol in that paper in order to consider limited 
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actuation and sensing update rates. A periodic event-triggered control technique for consensus 
of linear multi-agent systems is proposed in the present paper where each agent implements 
discretized and decoupled models of the states of neighbors. Communication delays are also 
considered. The case shown in lf44l represents only a particular case of the approach described 
here, when continuous actuation and continuous measurements are possible to implement and 
obtain, respectively. 

In the present paper decentralized event thresholds that guarantee practical consensu!] and 
strictly positive inter-event times are designed. The lower-bounds on the inter-event time in¬ 
tervals are independent of the particular system trajectories, therefore they hold for any two 
consecutive local events. The main contribution of this paper is the design of periodic event- 
triggered strategies for consensus of linear systems with limited sensing and actuation updates and 
with communication delays. The periodic event-triggered control technique automatically avoids 
the presence of Zeno behaviojj]. However, and for completeness, we establish the relationship 
between a selected sampling period and the performance of the consensus protocol with respect 
to the bounds on the state disagreement. 

The present paper addresses the several problems that were not considered by f[42l : 

1) Time-varying communication delays. The transmitted event-based measurement updates 
are subject to communication delays. 

2) Constrained sensing and event computation rate. Each agent does not need to continuously 
measure its own state but only at finite time instants. Similarly, continuous computation of 
errors and events is not necessary and these operations are performed only at discrete-time 
instants. 

3) Constrained actuation rate. This event-based approach also provides sampled actuation 
time instants instead of continuous actuation. 

Continuous measurement, actuation, and computation of events severely restrict the operation of 
the subsystems; therefore, the relaxations addressed in this paper offer a significant advantage 

'in the case of unstable trajectories practical consensus is obtained when the difference between the states of any two agents 
is bounded. 

2 Zeno behavior in event-triggered control refers to the occurrence of an infinite number of triggering events in a finite time 
interval. 
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in terms of implementation and resource management. 

The remainder of this paper is organized as follows. Section [Til provides a brief background 
on graph theory and describes the problem and the consensus protocol. Section [TTI] gives a result 
assuming continuous communication which will be used in the main results of this paper in 
Sections [TV] and [V] Design of periodic decentralized event thresholds for systems with limited 
sensing and actuation capabilities is addressed in Section [IV] Section [V] extends this approach 
in order to consider time-varying, but bounded, communication delays. Section [VT] presents 
illustrative examples and Section IVIII concludes the paper. 

II. Preliminaries 

A. Graph Theory 

Consider a graph Q = {V,£} consisting of a set of vertices or nodes V = {1,iV} and 
a set of edges £. An edge between nodes i and j is represented by the pair (i,j) e £. A 
graph Q is called undirected if (i, j) e £ (j,i) e £ and the nodes are called adjacent. 

The adjacency matrix A is defined by a tJ = 1 if the nodes i and j are adjacent and a tJ = 0 
otherwise. If (j, i ) G £, then j is said to be a neighbor of i. The set Mi is called the set of 
neighbors of node i, and N r is its cardinality. A node j is an element of M t if (j, i) E £. A path 
from node i to node j is a sequence of distinct nodes that starts at i and ends at j, such that 
every pair of consecutive nodes is adjacent. An undirected graph is connected if there is a path 
between every pair of distinct nodes. The Laplacian matrix £ of Q is defined as £ = V — A 
where V represents the degree matrix which is a diagonal matrix with entries da = V jeA! - a t] . 
For undirected graphs, £ is symmetric and positive semi-definite. £ has zero row sums and, 
therefore, zero is an eigenvalue of £ with associated eigenvector ljy (a vector with all its N 
entries equal to one), that is, £1^ = 0^. If an undirected graph is connected then £ has exactly 
one eigenvalue equal to zero and all its non-zero eigenvalues are positive; they can be set in 
increasing order Ai(£) < A 2 (£) < A 3 (£) < ... < Atv(£), with Ai(£) = 0. 

Lemma 1: Let £ be the symmetric Laplacian of an undirected and connected graph. Then, 
consensus is achieved if and only if 

1 / = x T tx = 0 , ( 1 ) 
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where £ 


£®Q,Qe R nXn is a symmetric positive definite matrix, 

r i T 


x(t) 


Xi(t) T X2 (t) T 


XN(t) T , 


Xi G R n , and ® denotes the Kronecker product. 
Proof. See Appendix. 


B. Problem Statement 

Consider a group of N agents with fixed communication graphs and fixed weights. Each 
agent’s dynamics are described by the following: 

Xi(t) = Axi(t) + Buftff), i = l,...,N, (2) 

with 

Ui{t fl ) = cF'^2{x i (t ll )-y j (t ft )), i = 1, ..., N, (3) 

jeMi 

where A £ R nxn , B £ R nxm , Xi £ W 1 is the state of agent i, and tq e M m is the control input 
for agent i. F e M’" xr ' and c £ R + are design parameters that are defined below. The variables 
Dj £ R n represent a model of the j th agent’s state using the decoupled and discretized dynamics: 

i) = Gyj(t M ), j = 1,iV. (4) 

for f /t £ where G = e A/l , /i = f /t+ , - and yj{t kj ) = Xj(t k .). We refer to /i as the 

discretization or sampling period since it is used to obtain the discrete-time model G and the 
state of each agent, x^ i = 1,..., N, is sampled every h time units. The discretization period h 
is constant but the communication intervals for each agent are not constant. The transmission 
time intervals are non-periodic and are determined by event-triggered rules. This means that 
at every time instant each agent samples its own state and updates its control input ©. It 
also uses its sampled state to compute its local state error and to determine if it is necessary or 
not to transmit the current state x(tf) to its neighbors. An example showing this relationship is 
shown in Fig. Q] This figure shows the continuous-time state xft) and its discrete-time model 
yftfj). Two event-triggered model updates are shown, at time t k . and at time 4 i+1 . The first 
event-triggered model update occurs at sampling time while the following event time instant 
is equivalent, in this example, to sampling time t M+ 7 , that is, 7 sampling time instants occur 
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Fig. 1. Representation of periodic sampling time instants and event-triggered communication time instants 


between and t^+i in this example. Note that at the event time instants the model is updated 
and it takes the current value of x % (t ^). 

The main advantage of using a periodic event-triggered control scheme is that measuring the 
state of each agent (and the associated computations that require evaluation of state errors and 
thresholds) is only performed at some periodic time instants (every h time units) instead of doing 
it continuously as it is common in most event-triggered control schemes. There are additional 
advantages of using discrete-time models of the continuous-time system dynamics compared 
to Il42ll . |[44l where continuous-time models were used. For instance, continuous actuation is 
not required and the operations related to implementing discrete-time models at every node are 
simplified with respect to implementing continuous-time models. 

Every agent in the network implements a model of itself y^t^) and also models of its neighbors 
Vj(t fj), j 6 Mi. The model state y,(f /( ) is not used by agent i for control since the real state, 
Xiftn), is locally available at every sampling time t^. However, the local model state, 
is used to trigger local events. This way of defining the local control input ([3]) represents a 
difference with respect to the approach in HTl and ll42l where the control input is a function of 
model variables, including y l . 

Local events for agent i are defined as follows. When an event is triggered, agent i will transmit 
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its current state Xi(t k .) to its neighbors and each one of the neighbors updates the model state 
they have of agent i. At the same time instant, agent i will update its local model using the 
same state measurement. The notation t ki represents the local broadcasting time instants, i.e. 
when agent i triggers an event and transmits its current state. Since we check for events not 
continuously but only at sampling times t then t k . is equal to some sampling instant but, in 
general, not every time instant t M is a broadcasting instant t ki . The decision to trigger an event 
in order to broadcast the current state is given by the event-triggered schemes described in the 
following sections. 

For negligible communication delays, agent i and its neighbors will all update their local 
models yi(t^) at the same time instant. Since agent i and its neighbors use the same measurements 
to update the models, say, Xi(t k J and the model dynamics © represent the decoupled dynamics 
where all agents use the same state matrix, then the model states implemented by agent 

i and by its neighbors are the same. In the presence of communication delays the previous 
statement will not hold and we will differentiate between ya(t M ), the model state of agent i as 
seen by agent i, and Urjit,,), the model state of agent i as seen by agents j, j E N r . More details 
concerning communication delays are presented in Section lYl 

The model update process is similar for all agents i = 1,..., N. The local control input © 
is decentralized since it only depends on local information, that is, on the sampled state of the 
local agent and on the discretized model states of its neighbors. Continuous or even periodic 
access to the states of neighbors is not needed. 

Note that the difference between ZOH periodic samples of the agent dynamics © and our 
proposed models © is given by the input term in © and this input decreases as the agents 
approach a consensus state. It can also be seen that in the particular case when systems © 
represent single integrator dynamics, then our models degenerate to ZOH models as in IfTTIl . 

m. 


III. Consensus with continuous measurements 

This section provides an important result in Lemma [2] that will be useful in the following 
sections. Let us assume in this section that continuous actuation by each agent and continuous 
communication between agents are possible, then the agent dynamics © are now driven by the 
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local control inputs 


Uj(t) = cF (xj(t) - Xj(t)), i = l...N. ( 5 ) 

jeAfi 

Assume that the pair (A, B) is controllable. Then, for a > 0 there exists a symmetric and positive 
definite solution P to 

PA + A t P - 2 PBB t P + 2 aP < 0. (6) 

Let 

F = -B t P (7) 

and c > 1/A 2 . 

Theorem 1: Assume the pair (A, B) is controllable and the communication graph is connected 
and undirected. Define F as in ([7]) and c > 1/A 2 . Then the following symmetric matrix 

C = CA c + A T c C (8) 

has only n eigenvalues equal to zero and the rest of its eigenvalues are negative. In addition, 
the eigenvectors associated with its n zero eigenvalues belong to the subspace spanned by the 
eigenvectors associated with the n zero eigenvalues of £, where £ = £ ® P, A c = A + B, 
A = In <£> A, B = cC® BF, and 7y is an identity matrix of size N. 

Proof. See Appendix. 

Lemma 2: Assume the pair (A, B ) is controllable and the communication graph is connected 
and undirected. Then, protocol ©, with F defined in ([7]) and c > 1/A 2 , solves the consensus 
problem for agents described by ([2]). Furthermore, the Lyapunov function defined by V = x' C.i: 
has a time derivative along the trajectories of © with inputs © given by V — x T Lx. 

From Theorem Q] it can be seen that V is negative when the overall system is in disagreement 
and is equal to zero only when the corresponding states are in total agreement. In the latter case 
we also have V = 0, see Lemma Q] Different from consensus with single integrators, where 
the agents converge to a constant value, here it is only required that the difference between 
states of agents tends to zero, regardless of the particular response of the systems. As with many 
consensus algorithms, an estimate of the second smallest eigenvalue of the Laplacian matrix is 
required; this is the only global information needed by the agents. Algorithms for distributed 
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estimation of the second eigenvalue of the Laplacian have been presented in Ii45ll . If46l . Readers 
are referred to these papers for details. 


IV. Decentralized Periodic Event Triggered Consensus 

In this section we consider the case when agents use event-triggered communication strategies 
in order to reduce the frequency of transmissions. Every agent implements discrete-time models 
© and samples its own state every h time units in order to compute its local control input and 
check its local thresholds. It is assumed in this section that communication delays are negligible. 
We derive decentralized thresholds that depend only on local information and can be measured 
and applied in a decentralized way. Additionally, it is shown how the choice of the period h 
affects the bounds on the disagreement between states of agents and how to select the threshold 
parameters to determine lower bounds on the inter-event time intervals. 

Define 


e*(t) = Vi(t) ~ Xi(t ) 


(9) 


and x = 


T 

x\ 


T 
• X N 


y = 


T 

yi 


T 

■Vn 


e = 


T T 

e 1 ... e N 


, where yi(t) represents the hypothet¬ 


ical continuous-time model 


Vi(t) = Ayiit) i = 1...N. 


( 10 ) 


Definition (fTOl) is only used for analysis of the event-triggered controller; it is not used for 
implementation of the event-triggered consensus algorithm. 

Let us define the discretization errors 

Xi ( t ) = Xiit ^) - Xi ( t ) 

Vi ( t ) = Vi ( tp ) - yi ( t ) (11) 

e;(t) = eiit^) - a{t) 


for t E [tf^tn+i). Also define x = 

vj 1 oj 1 

. . . Xjy 

T 

> y = 

yi ■■■ Vn 

T 

, e = 

e{ ■■■ e N 

c, the new coupling factor c\ = 2c is now used in the inputs d3 

). 



Theorem 2: Assume the pair (A, B ) is controllable and the communication graph is 
and undirected. Define F in 0 and c\ = 2c where c > I/A 2 . Then agents © with 


Instead of 

connected 
inputs © 
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and discrete-time models © achieve a bounded consensus error where the difference between 
any two states is bounded by 


lim^oo \\ Xi {t) - Xj(t)\\ 2 < p ^ P) 

for i. j = 1 ,N and for any t > 0, if the local events are triggered when 

5 i > ac\zjPBB T Pzi + 77 , 


where 0 < a < 1 , /3 = 


^min^O ( X) 


Amax(^) 


> 0 , 


b > 0, 


and 


<5, = cN^N - 1) + ^ )ef(jt /l )PBB T Pe i (t li ) 
+ ci(l + 2bNi)zf ( t~)PBB T Pzi(t~) 

+ cN z (^ + 3b(N-l))er(t-)PBB T Pe l (t-) 


* = ( X A-vAh)) 

jeNi 


Zi = & - vA 

j&K 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


Furthermore, the agents do not exhibit Zeno behavior and the inter-event times t^+i — tki for 
every agent i = 1,N are bounded as follows 


0 /i <C + l i/,' ( 


(17) 


if 


r 1 > (2c 1 7V i (6(AT - 1) + f) \\E|| 2 + Cl (l + 2bN l )(2 + b e \\E \\) 2 ) \\PBB T P\\ z, 2 (18) 

where 

Zi = A max (£)w Vm , (19) 

V A min (£) 

V M = max |l/( 0 ), ^}, E = f Q h e A ^cBFds, b e = ^N,N(N - 1 )(| + i), and h = t^-t^ 

is the discretization period. 
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Proof. By implementing the coupling factor c\ = 2c in 0), we can write ©-© in compact 
form as follows: 


x= Ax + B v x(tfj) + B A y(t M ) = (A c + B)x + B A e(t + B x x 


( 20 ) 


where B v = c{D®BF , fE = — Ci*4.(g)-BF, i?! = c x C®BF. The overall system dynamics (l20l ) 
is written in terms of the closed-loop state x plus two error variables that are introduced in the 
system dynamics because of the sampled inputs and the event-based communication strategies. 
The error eftf) is due to the fact that each agent only communicates its local state at some local 
event time instants and the error x is due to the fact that only discrete models of neighbors and 
periodic samples of the local state are used to compute the control input instead of continuous 
variables. In other words, the first error results from limiting the communication between agents 
and the second error results from reducing actuation time instants at each local node. 

Consider the candidate Lyapunov function V = x T Cx and evaluate the derivative along the 
trajectories of systems © with inputs ©. We can express V as follows: 

V = x 1 C {(A c + B)x + B A e{tfj) + B x x) + {{A c + B)x + B A e(t lt ) + B x x) T Cx 
= x 1 Cx + 2 x T CBx + 2x J CB A e{tfj) + 2x T CB x x 

= x T Cx + 2 Eli EfceA/ifo - x k ) T PBB T P (21) 

(- c'EjetfM ~ x i) + Cl EjeAfi e j(^) 


-ciE 


jeMi 


[Xi - Xi 


Eq. (HIT) can be further written in terms of the sampled state information xft^), Xjitf) and of 
the discretization errors Xj defined in (fTTl) as follows: 

y= xT Cx + 2Ei, E“ x k)) T 

PBB T p(-c'£ jeJsf .(x i (t li )~ x j (t li ) - ( Xi-Xj )) 


c i JCjeMi e j(^r) c i Ej 


jeJSi 


[Xi 


Xn 


= x T Cx + 2Ef =1 [E 

PBB T P( -cE ( x * fa) ~ x i (tn)) 


( 22 ) 


where 


+ c i YljeM e A t r )) + 6 


Hi=T. k ^.{ii-iil T PBB- r P[c 


X P 


' C 1 Yhj^ATi e j(t, 


(23) 
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By using the state error definition © and the definition of z r in (fT5l) we can write (l22l) in the 
following form: 


V = x T Cx + 2 Eir [(* + ^(t,)) T 

PBB T P^-c(Zi + e j(tn))+Cl Y^jeAfi e 3^S) + & 
= x T Cx + 2 El, - czjPBB T Pzi 

+ e T k (t,)PBB T PZ j€ M + & 


( 24 ) 


Using the inequality | x' y | < + j^y 1 y, for 6 > 0, we have that the following expression 

involving the second term in (1241) holds 


EfceM e K^) PBBTp EjeVi e A 


< 


Y.mApb&pY., 


jeJVi e J 


<Y,^,T,i<MM PBBTp 4\ 

<N4 + i) E i ^,eJ p BB Tp e j 

Since the communication graph is undirected we have the following properties 

TV TV 

EE eJ PBB T Pej = EE eJ PBB T Pei 

i= i jeM i=i jeM 

and 


(25) 


(26) 


Eli JVi EjgM eJPBBTPe, 

= Ell E EjeVi e fPBB T Pei (27) 

— Ell(^ _ 1) EjGTV; e J PBBTp>e i 

Then, using (l27l) . we can write the following expression which bounds the second term in (l24l) 


Eli N i(\ + Y b ) EjeM eJ pB B Tp ej 
= Eli E(! + i) E^ eJPBB T P ei ( 28 ) 

< Eli - 1) (| + i) efPBB T P ei . 

Let us now analyze the term (1231) . Using the definitions (fTTl) and (fl6l) the summation of terms 
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& can be written as follows 


czj PBB t Pz; 


E N £ _ y^iV 

i=1 l^i=1 

+ci z T i PBB T PY J 


jeMi e i 




jeM e i 


(29) 


- Cl zJPBB T PY, 3&K e 3 {t^ 

+C 1 Y. k ^ PBBT P^MM^ 


We again use the inequality ary < |ar x + ^jry, for b > 0, to obtain: 

Eii6< Eil [c(l + 2 bN i )zJPBB T Pz i 
+ cNi (^ ^ 


Using 


+ 2^zi2)gT PSS Tpg. 

and (QOl) we can bound ([24b as follows: 


(30) 


U< x T £x + Eli 


c P , 

3W-1 


T PBB T P Zi 


cN t (b(N- 1) + ^t±)eJ(t,)PBB T P ei {t,) 


(31) 


+ ci(l + 2bNi)zf PBB t Pzi 

+ cN i (^+3b(N~l))%PBB T Pe i 

Note that we can compute errors and make communication decisions only at the sampling time 
instants that is, only when each agent measures its local state. Also note that the variables 
Zi and e-i are equal to zero at every time regardless if an event is triggered or not, i.e. 
z^t/j) = ej(f M ) = 0, for /i — 0,1, 2,.... However, it is necessary to consider the non-zero value 
(before they are reset to zero) of these errors in order to evaluate ([3Tb and to decide if an event 
needs to be triggered or not. In other words, at every time instant we need to evaluate event 
thresholds using the discretization errors just before they are reset to zero, that is, just before 
the actuation update occurs at time t 

Let us denote the time just before the actuation update as t~, so we check the discretization 
errors in (OTT) using Zi(t~) and ej(t“), i.e. the errors just before they are reset to zero. Thus, we 
can write 


N 

V < x T Cx + [ - c 1 zfPBB T Pz i + <SJ 

i= 1 


(32) 
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Time 


Fig. 2. Computation of discretization errors based on current and previous measurements 


where - &(*„)), and the variables x^), 

and Vj{t~), for j e M, can be obtained locally (see Fig. [2]) at every time using the 
current and previous samples of the local state and the local models as follows 


Xi(t^) = Xi(t^- 1) - Xiitfj) 

= 3/jM 


(33) 


The variables Xj and yj, j = i,j E Afi, cannot be measured continuously but only at times f /7 as 
it is illustrated in Fig. O 

When threshold (fl3l) holds, say, at time t ki then the error resets to zero, that is, ei(t ki ) = 0, since 
Uiitki) = Xi{t ki ). Also, since t ki = t,, for some non-negative integer /i (recall that communication 
events can only occur at the sampling times t^), then, the control input is updated and the 
discretization errors are equal to zero, that is, z^t^) = e^t^) = 0. This means that <),(/:/,J = 0, 
then the following holds 


Si < ac-izf PBB t Pzi + 7], 


(34) 


Consequently, 


N 



(35) 
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Since a — 1 < 0 we have that (a — l)cizf PBB T P Zi < 0 for % — 1, N and 

V < x t jCx + Nr] (36) 

We use the fact that C is positive semi-definite and it has n zero eigenvalues with corresponding 
eigenvectors v±... v n . Let x = x\ + x 2 such that (xj,x 2 ) = 0 and Cxi = 0, that is, x\ belongs 
to the subspace spanned by v 1 ... v n . Consider 

x t jCx = ( x\ + X 2 Y£(x\ + X 2 ) = X 2 CX 2 (37) 


also 


x[Lx2 < X max (C) X 2 X 2- (38) 

From TheoremlHwe know that — C is positive semi-definite with n zero eigenvalues and Cx\ = 0. 
We can see that 


and 


x 1 Cx = ( x\ + X 2) 1 C{x 1 + X 2 ) = X 2 CX 2 


(39) 


xl(~C)x 2 > XrmnM-£) X 2 X 2 


(40) 


where A min ^o(—£) > 0 is the minimum eigenvalue of — C other than zero. Combining the 
expressions above we have 


X2CX2 < 


^min^OV T n 

-X 2 CX2 • 


Amax(-^) 


( 41 ) 


Then V can be bounded as follows: 


V < -Px T Cx + Nr] = -PV + Nr]. (42) 

Solving (1421) we have that 

V{t) < e-^C(O) + Nr] f* e-M^dr 

< (V( 0 ) - e~' 3t + (43) 

Expression (1431) represents a bound on the consensus states as a function of the initial separation 
of the agents V(0) = x(0) r £x(0). We can express V(t) — \ Yln=i J2jeJsS Xi ~ x jYP{ x i ~ x j ) 
and a direct bound on the difference between any two states i,j can be obtained as follows. 
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Since the graph is undirected the term (x t — Xj) T P(xi — Xj) appears twice in the summation 
V(t), then we can write 


Amin(-P) II Xi - Xj || 2 < (Xi - Xj) T P(Xi ~ Xj) 

) \ (44) 

Finally, the difference between any two states can be bounded as in (fl2l) . 

We will now prove that the inter-event times are lower bounded by the sampling period as in 


First note that 


where C n 
holds 


*ill = 

E jeM Pi(p) -^(p) - 



< 


+ 

YhjoNi 


— Il'^n^'ll ll* / ^-^^|| 

£ <g) I n and A n = A 0 I n . For the first term in (1451) we have that the following 


I £71^11 y X~x P n P n x ^ A max (C n )Vx T X 

c (£). 


< A 

< A 


_ / 'max V 


-^min(-£) 

Vm ~ 

^min(£) 


Te second term in (1451) can be written as follows 


(46) 


||*^-n^|| yfe AnA n C 

= \JtA± Slew e l 

< yjzli mr - 1 )(! + i)efe, 
— b e || ('i || 


Then we can write (1451) as follows 


(47) 


pill < A max (£)y^AAi_ + b e IpH (48) 

for any time p. In particular, at the local event times t k . we have that e r {t ki ) = 0 and 11 z, t ( t k .) 11 < 

Zi- 

In order to prove (fT71) we need to guarantee that (IT4l) does not grow from zero at time t ki to 
r] at time t k . + h and will not trigger an event at t k . + h. In other words, for a given sampling 
period h, we need to find the value of r] that guarantees (fTTT) . 
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Consider the error dynamics of agent i from t ki to t ki + h as follows: 


ei(t k . + h) = Ui{t ki + h) - Xi(t ki + h ) 

= Gyi(t k .) - Gxi(t ki ) - Ezi(t ki ) (49) 

Gei(t ki ) E Zi(t k 2). 

Note that ei(t ki ) = 0 since the error is reset at the event time t ki , then we have that 


ei(t ki + h) = —Ezi(t ki ). (50) 

The effect of the discretization errors is considered before they are reset to zero at time t ki + h 
because of the actuation update, so we have 


e-i{t ki + h) = ei(t ki ) - ei(t ki + h) 

= —e-i(t ki + h) = Ezi(t ki ). 

Now, let us analyze the term Z{. Consider 

P i{tki + h)\\ = ^2 jeM .(xi(t ki + h)- Xj(t ki + h)~ e.j(t ki + h)) 
< || C n x(t ki + h)\\ + b e ||e i(t ki + h )|| 

Note that 


(51) 


(52) 


\E n x{t ki + h )|| = || C n x(t ki ) - C n x(t ki + h)\\ 


< 


< 


-^%-AVvK) + Vv(th + h)) 


(53) 


2A n 


c(C) 


yj Amin (C) 


A JVm 


We use ®, and (l53l) to analyze the growth of the term 5 t from t ki to t ki + h 


Si(t ki + h) < 2ciNi(b(N — l) + f) \\E || 2 \\PBB t P\\ z?(t ki ) 

+ Cl (l + 2WVj) ||P55 t P|| ( ^ II^H ^(4j ) 2 

\VA m in(>C) / 

Then, by the selection of r/ in ([J~ 8 l) we can guarantee that after an event instant at t k . the term 
Si cannot grow from zero to 77 and threshold (fl3l) is not triggered. Then, the inter-event times 
are greater than h > 0 and, obviously, we can guarantee that Zeno behavior does not occur at 
any node. • 

Remark. The design parameter 77 provides a tradeoff between performance as measured by 
the consensus error (fl2l) and the bound on the inter-event intervals (fT71) . Also note that the 
variables used to compute the threshold (fl3l) . which define the events at node i, are available 
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locally. Concerning global information we only need an estimate of the second eigenvalue of 
the Laplacian, as it was mentioned earlier. 

Remark. Note that the triggering of a local event by agent i, i.e. the transmission of a 
measurement Xi(t k .), does not change the local control input u, beyond the local actuation 
update, since v n is not a function of y t . The transmission of a new measurement updates the 
control inputs of neighbor agents and for the local agent, it only resets its local state error t\ 
and Si. 

Remark. The results provided in this paper hold for agents described by general linear dynam¬ 
ics. The common cases of single integrators and double integrators are particular cases covered 
by this framework. The single integrator is modeled as a ZOH and the double integrator is 
modeled similar to lfl9l . that is, velocity as a ZOH and position as a first-order-hold model. 

Remark. The event-triggered consensus algorithm provides a bound on the difference between 
any two states not only at the sampling time instants but for any t > 0. Also note that rj can be 
made arbitrarily small (thus decreasing the bound ( [Til) ) and still satisfy (fT71) by choosing a small 
discretization period h since we have that both E —> 0 and I — G—t 0 (therefore S t {t ki +h) —> 0) 
as h —> 0. 

V. Decentralized Periodic Event Triggered Consensus With Communication 

Delays 

In this section we consider the presence of time-varying but bounded communication delays. 
Since the measurement updates will be delayed, the neighbors of an agent i will have a version 
of agent V s model state that is different than agent V s version. It is necessary to distinguish 
between the model state as seen by the local agent, itself, and as seen by its neighbors. Define 
the dynamics and update law of the model state of agent i as seen by agent i as 

yu{tfi+1 ) Gy a (^)) yn(tki) Xi(t kj ^. (55) 

The measurement Xi(t ki ) is transmitted by agent i at time t ki and will arrive at agents j, j E A/), 
at time t ki + di(t ki ). For a given update instant all receiving agents experience the same delay 
di{t k .). However, this is not a constraint and the communication delays can be generalized so, 
for the same update instant t k ., the neighbors are updated at different time instants. This is 
described at the end of this section. 
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Assume without loss of generality that di(tp) = hpi(tp) and p % > 1 is an integer, that is, the 
delay is an integer multiple of the sampling period h. If the delay is not an integer multiple of 
the sampling period the receiving agent uses the corresponding delayed measurement at the next 
sampling period to update the corresponding model which effectively makes the delay to be an 
integer multiple of h. 

Define the dynamics of the model state of agent i as seen by agent j, j e J\f t , as 

yij(tfi+i) Gy % j (t/j) i (56) 

The states of these models are updated when a delayed measurement of agent i is received by 
agent j, j e J\f % . The state measurement x % (tp ) transmitted by agent i at time instant tp is 
received by agent j, j e Hu at time instant tp + hpi(tp). Let us define the update law of the 
model state of agent i as seen by agent j, j e J\fi, as 

ytjitp + hpi (tp )) = fd(xi(tp),pi(tp)) (57) 

where tp represents the update instants triggered by agent i and hp,{t k .) represents the commu¬ 
nication delay associated to the triggering instant tp. 

Define a positive and constant upper bound on the communication delays by d = ph < 
tp + 1 — tp, that is, di(tp) < d for any triggering instant tp and for i — 1, N. Later in this 
section we will define the design parameters that bound the inter-event times as a function of 
the delay d. Note that p > Pi(tp), for i = Assume that the current delay, Pi(tp), is 

known to the receiving agents, for instance, by applying time-stamping techniques. 

Since both, yu and y l; j , use the same state matrix to compute their response between their 
corresponding update instants, then we define 

fd(xi(tp), hpiitp)) = Xiitp) (58) 

that is, the delayed measurement is propagated forward in time and the result is used to update the 
state of the model as shown in (l57l) . By definition, we have that the following local triggering 
event will occur at time £ fc . +1 > tp + d, this means that y tl {tp + di(t k J) = G p,l ' t p'>x t (tp) 
because no other local event has been triggered since time instant tp. Therefore, we have that 

yaitf,) 7 ^ Vijitfj) for e [tp,tp + hpi(tp )) and yuit^) = y^t M ) for t e [tp + hpi(tp),tp +1 ). 
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Fig. 3. Relation between state Xi, model states ya, ytj, and corresponding errors ea, e (J 


Define the state errors 


Unity ,) 

^iji^y) Uiji^y) •^ii^y)- 


(59) 

(60) 


Note that eu{t ki ) = 0 and e*^) = e,,^), for e [t ki + hpi(t ki ), t ki+ i). These relations are 
pictured in Fig. 0 


Also define C (ty) = e^) ... e^ JV (f M ) 


1 T 


r>nN 


and ( d (ty) — ef jl (ty)... e Nj N (ty) 


where the components e tJ> in Q represent the errors defined in (l60l) . that is, the error of 
agent i as seen by its neighbors j l , ji e TV). 

Define the discretization errors y lt [i) = yu(t^ - yu(t), y l3 (t) = y^t M ) - y t] {i), e vi {t) = 


e ii{ty) e ii(f)j e ij{ty) V Ui1---Un1 V 


1 T 


, Vd = 


yln---y T N jN 


, c = 


'-'T 

e ll ■ ■ ■ e NN 


> C d 


Kj r r 

e lji ■ ■ ■ e Nj N 


The dynamics of every agent in ([2]) with communication delays captured by the new input 
definitions (which are functions of delayed model states yp,(ty)) : 


Ui(ty) = CxF (Xi{ty) - Vji(ty)), % = 1...N, 
j&Mi 


( 61 ) 
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can be written in compact form as follows: 


x = (A c + B)x + B A Q(t^) + B ix 


(62) 


where the coupling strength c\ = 2c has been used. 

Note that if we follow the same analysis as in Theorem [2] we can arrive at an expression similar 
to (f3J~b involving the local control inputs z r , the local discretization errors, and the delayed errors 
e tJ (in ). The latter will create a major difficulty in designing the local events since the local agent 
i does not have access to the errors as seen by its neighbors. More importantly, the local 

agent is not able to reset the error eij(t but only the local error eu(t^). 

The following theorem provides a method to design local events in the presence of commu¬ 
nication delays and using the local state errors to locally evaluate the events. 

Theorem 3: Assume the pair (A, B) is controllable and the communication graph is connected 
and undirected. Define F in © and c\ = 2c where c > 1 /A 2 . Then agents ([2]) with inputs (16TI) 
achieve, in the presence of communication delays di < d, a bounded consensus error where 
the difference between any two states is bounded by (fl2l) for i , j = 1,.... A r , if the events are 
triggered when 


Si > (TCizJPBB t Pzi + r] 


(63) 


where 0 < a < 1, 


Si = Ci (1 + b) 2 zJ{t^)PBB T Pzi(t M ) + 5 di + Si. 


(64) 



(65) 


jeMi 



Si = Cl (l + 6)(1 + i)A|e£(g((G - I)G P ~ 1 ) T (G - 
+ 2 \\(G—I)G p ~ l ea (f M )|| (||G-/|| T h +\\E\\)zi 
+ ((\\G-I\\T h +\\E\\)zi) 2 }, 


( 68 ) 
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where b> 0, A = \ max (A 2 ® PBB T P), z t = &e = /^(^ “ !)(| + ^)> 

T = / Q d || e A ( d - a )c5F|| cis < £, T h = f*~ h \\eA d - h - s ) c BF\\ds, E = f Q h e A ( h ~^cBFds. 

Furthermore, the agents do not exhibit Zeno behavior and the inter-event times t ki +i — t ki for 
every agent i — 1, N are bounded by d > 0, that is 

0 < d < t ki+1 - t ki . (69) 


if 


(70) 


, > Cj(l + bf \\PBBT P \\ + TO 

+ e 1 (l + i)(||G»|| + l) 2 T 2 5? 

+ c 1 (l+6)(l + i)(||(G-/)G»- 1 ||T + ||(G-/)||T k +||ii;||) 2 2 f 

Proof. Consider the candidate Lyapunov function V = x 7 Cx and evaluate the derivative along 
the trajectories of systems © with inputs (f6TI) . 

V = x T C (( A c + B)x + B A C d [tr) + ^x) 

+ {{A c + B)x + B A Cd{tii) + B\x) Cx. 

Following similar steps to (l2TT) - (l24l) in the proof of Theorem [2] and using the inequality |k T 2/|| < 


(71) 


bJT n 


x 1 x + j^y 1 y, for b > 0, we can write 


V<x T Cx + 


- czf PBB T Pzi 


c (l~f b ) e ki 

+ c(l+6) T. kM (ii-i*) T PBB T PY. i 


(72) 


jeAfi 


-Xj) 


where ejft^) represents the state error, at time instants t,,, of agent j as seen by agent i, for 
j e M;. We can write (l72l) as follows 


1/ < x T Cx + ffZi I - c 1 z[PBB T Pz i 

+ Ci(l + i) Ysk&K ^2jeATi e ji{ t r) 

+ ci (1 + 6) (zfPBB T Pzi + 2 z?PBB t PJ2 


(73) 


jeAfi e P 


E k ^ i e T ki PBB T PE je M i e j i 


23 













Using the inequality ||o; T ?/|| < \x T x + h-y T y, for b > 0, once again, we obtain 


2 ^ ^ 1 2b 
T r — i 


V < x T Cx + - c l zfPBB T Pz i 

+ ci(l + b) 2 zJPBB T Pzi 

+ ci(l + |) e ki(^)P BBT P Eje Mi e i^%) 

+ c 1 (l + b)(l + l)E k ^ i e T ki PBB T PT lJ 


jeM e i i 


= x T Cx + Ef=i 


c\zj PBB T Pzi 


+ ci(l +b) 2 z[PBB 1 Pz, 


(74) 


+ c 1 (l + l)C(t ll )(A®B T P) T (A®B T P)C(t fJ ,) 

+ Cl (l+b)(l+l)(J(A®B T PnA®BTp)Cl 

The variables Zi and z r can be computed locally by every node. Let us then focus on the delayed 
error terms Cj(^) and Cj • Define v^) = y^t M ) - yn(t^, then we have that 


&ij (f/i) yijit/j) 

Vij(h) (.yuit/ji) (75) 


For the term ty(f /t ) the following holds: 




Gu i (t lu ), tp e [tkiAki + di(t k ..)) 

0, G [tfej + ^i(ffei)) ffcj+l) 


(76) 


with Vi(t k .) = e it (t k ). The notation t k , represents the event time t ki but just before the local 
error is reset to zero. This update of the variable v at time t ki is obtained by simply realizing 
that Vi(t ki ) = yij(t k .) - ydt kt ) = y n (t kj ) - x t (t k .) = e tl (t kt ) (the local error just before it resets 
to zero), since the local model y tl is updated using x r {t ki ) and yij(t ki ) = yij(t k .) = yu(t ki ). 

Define = [uf (f M )... and we have that 


CaM = K^) + C(g- (77) 

Consider the worst case scenario (greatest difference between ( d (^) and ((t^)) given when all 
agents transmit at the same instant t k and the greatest possible delay d (> hj>, ) is present. Then, 
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we have the following 


Cj(4 + d)(d(tk + d) — ( T (t k )(G p ) T G p ((t k ) 

+ 2C T {t k + d)G p ((t k ) + C T (4 + d)((t k + d) 

= + 2e l(^i + d)G p en{t k ,) 

+ e !(^i + d)eu(tki + d)] 

- Sill [ e S(4 i )( ( ^ P ) T ^ Pe **( i fe i ) + 2 || GPe ii( i fc i )|| + (T^, m ) 2 ] 

(78) 

where G = I N ® G. The local error efft'y 7.) represents the error just before the update instant, 
i.e., before it is reset to zero because of the update at time t ki . On the other hand, e lt (t ki + d) 
represents the error after the update at time t k . and it can only be estimated using 

\\eu(t ki +d)\\ < ||e Ad eii(f fei )|| + z i>M Jq \\e A< ^ d ~ s) cBF\\ ds 

< 'Tzi,M 

where ea{t k .) = 0 and z h M represents an upper-bound on the the local control input for the time 
interval e [t ki ,t ki + d], that is, || Zi(t ki + hpi )|| < z i>M for = 1, ...,p. 

Since the worst case is given by the maximum delay d we can use (1781) and the current local 
error eu{t^) to bound the delayed error + d) for any sampling time ^ > 0, therefore the 
term (l67l) is used as a part of the overall threshold (l64l ) . In other words, we propagate the current 
error eafyf) forward in time using the worst case delay d, as if we had an event at time t ll: , then 
we check (l63l) . If (1631) holds then an event is triggered; otherwise, no event is needed and we 
repeat the same process at the following sampling time instant 

Now, we consider the discretization error Q corresponding to the delayed state errors ( 

The discretization errors are reset at every sampling time t therefore, we only need to consider 
the effect of these errors at time t k + d for the worst case delay d (> hp, ). In other words, we 
need to evaluate the difference (d(t k + d) = (d(t k + d — h) — Cdifk + d) as a function of the state 
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errors ((t k ). This can be done as follows: 

C d{tk + d) T Cd(tk + d) 

— ((d(tk + d~ h) — (d(tk + d)) T ((d(tk + d — h) — (d(tk + d)) 

= c(xn& - - G^)at-) 

+ 2 (((tk+d) — ((t k +d— h)) T (G p — G p X )C (t k ) 

+ (((tk + d) — ((tk + d — h)) T (£(tk + d) — ((tk + d—h)) 

= El 1 [&kM G - I)G p ~ 1 ) T (G - I)G p -'eu(t- ki ) (80) 

+ 2 (ea(tki+d) — eu(t ki +d—h)) T (G — I)G P 1 eu(t k ..) 

+ (&u(tki+d) — en(t k .+ d — h)) T (en(t k .+ d) — ea(t ki +d—h))~\ 

= Eli [4(^)((G - I)G p ~Y(G - I)G p ~ l ea(t kj ) 

+ 2 ||(G-/)G^(4-)|| (||G-/|| T.+ PUKm 

+ ((l|G'-/||T ft +||£;||)^ m ) 2 ’ 

since we want to guarantee (l80l) for any sampling time then we use the current error e^(f M ) 
and we include (1681) in the overall threshold 
Then, the following holds: 


N 


V < x T Lx + [-ciz?PBB T Pzi + <Ji] . 


(81) 


i= 1 


Then, the local thresholds can be defined based on the local errors ea(t^) as in (l63l) with d; 


(82) 


given by (1641) . When an event is triggered the error e# is reset to zero and the following holds 

V < x T Lx + ci Eii(^ - 1 )zJPBB T P Zi + Nr] 

< x T Lx + Nr] 

and the bound (fl2l) on the difference between any two states follows. 

Also note that 

|| Zi(t ki + hpi )II = YjjeK ( x i(tki + hpi) -Xj(t ki + hpi) - + hpi)) || 

^ || ^"nX{tki ^Pi)|| H - be \\^ij (j-'ki H - ^Pz)|| 

< A max (£)+ b e || Zi (t ki + h Pi )|| ft* \\e A ^-^cBF\\ ds 

< ^(L)^^^ + b e || z^ + h Pi )|| f 0 hpi \\e A ^-^cBF\\ ds 

Expression (l83l) can also be written as 


(83) 




A m in(-£) 


(84) 
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From (1431) we have that max ((i { Vit^} = V M = max |v(0), ^ j and since the right hand side 
of (l84l) increases as p, increases from 1 to p, then we can conclude that 


Zi{t ki + hpi) || < Zj, tf = Zi 


(85) 


for pi = 1, ...,p. 

The final task is to determine r/ such that the inter-event times satisfy (1691) for a given d that 
satisifies T <b At time t k . we have that the local error is reset, i.e. e lt {t ki ) = 0, and no event 
is to be triggered by the local agent i during the time interval t /( e [t ki , t ki + d). Then, the error 
ea(tki + d) can be estimated as in (f79l) . Now, the term S l {t, ki + d) can be bounded as follows: 


\\S(t ki + d)\\ < Cl (l + b) 2 zf{t ki +d)PBB T Pz i (t ki +d) 

+ Cl (l + i)(||G^| + l) 2 T 2 ff (86) 

+ Cl (l+6)(l + i) ( ||(G—J)G P_1 || T + ||(G-/)|| T h +\\E\\ )'zf 
Following similar steps as in (l83l) the term Zi(t ki -\-d) can be estimated as follows: 


\zi(t ki +d)\\ = £ ie A/i + -Xj(t ki +d) - eji(t ki + d))\ 

< || C n x(t ki + c0|| + b e ||e»j(4i + d) II 

< A max (£)y^ ^ m .^ ) (rrTf + i-b e T fc ) 


(87) 


Then, by the selection of r/ in (f7Ql) we can guarantee that after an event instant at t k . the term 
Si does not grow enough to reach the value r/ at time t k . + d and threshold (1631) is not triggered 
before or at t ki +d. Then, the inter-event times are greater than d as shown in (l69l) and, obviously, 
we can guarantee that Zeno behavior does not occur at any node. • 

The communication delays are time-varying, that is, the value of the delays Pi(t ki ) and 
Pi(t ki+ 1 ), for any are not the same in general. However, it was assumed in this section 
that for a given update time instant t ki the neighbors of agent i will all receive the measurement 
Xi(t ki ) at time t ki + hpi(t ki ). This can be generalized to consider different delays from agent i to 
each one of its neighbors at the same update time instant. In other words, we can consider the 
delays Pij(t ki ) and Pij>(t ki ) to be different in general, for j,j' e A/). Thus, the models y tJ and 
'ijiji may differ during some time intervals, when one is updated sooner than the other. Fig. [4] 
shows a simple example where agents 1, 2, and 3 are neighbors of agent i. In general, we have 
that pn(t ki ) ^ Pii{t ki ) 3 Pii (tk,, ) and the delayed models are updated at different time instants. 
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v^p.i y- h p,2 v +h p i3 


Time 


Fig. 4. Relations between model states and corresponding errors when different communication delays Pij{tki) occur at the 
same update instant 


Assuming that p t j (t ki ) < p for j e M,, then it is clear that the same results in this section apply 
to the more general case under discussion since the growth on each error is still bounded by 
the error difference propagated using the upper-bound on the delay d = hp > hpij(t k ,), for any 

j e Mi. 


VI. Example 


Consider a decentralized model-based implementation of four second-order agents (N = 
4, n — 2) with unstable linear dynamics given by: 


Solving © we obtain 



0.2 -0.8 


0.7 

A = 

0.26 0.05 

,B = 

-1.1 


P = 


0.5859 -0.1575 

-0.1575 0.4274 


The continuous-time state matrix is unstable with complex eigenvalues A(A) = 0.125 ± .5 i. The 
discretization period is h — 0.002 seconds, the delay bound is d — 0.014 seconds, and we 
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compute r] = 10.85. The nonzero elements of the undirected adjacency matrix are a\ 2 = a 2 3 = 
a 34 = 1, (the corresponding symmetric elements are also equal to one). The initial condition of 
each agent are as follows 



-5.5 


1 

1 

07 

£ 1 ( 0 ) = 

- 6.1 

,z 2 ( 0 ) = 

-1.5 



5.9 


12.35 

£ 3 ( 0 ) = 

2.5 

,x 4 ( 0 ) = 

15.1 


Fig. [5] shows the response of the agents where it can be seen that the agents synchronize their 
states in each one of their two dimensions. Communication delays for each agent are time- 
varying and they take random values form the finite set: {0.010, 0.012, 0.014} seconds. Fig [ 6 ] 
shows the time instants where each agent broadcasts a measurement. Here, it can be seen that 
agents transmit information less frequently as they transition into a consensus state. The inter¬ 
event times for every agent are lower-bounded by the delay bound, that is, t ki +i — t ki > 0.014 
seconds. Due to a higher density of transmissions at the beginning of the simulation, a new 
figure, Fig. [71 has been included to show with more clarity the first 3 seconds of the simulation 
and to show that the inter-event time intervals for any agent are always greater than the delay 
bound d. 


VII. Conclusions 

Synchronization of state trajectories of linear multi-agent systems was studied in this paper. 
Multiple issues affecting the convergence to common trajectories were considered such as: limited 
sensing and actuation capabilities, limited communication, and time-varying communication 
delays. Event-triggered control schemes were proposed in this paper which not only provide 
decentralized control inputs but also allow for decentralized design of transmission instants where 
each agent decides, based only on local information, when to broadcast its current measurements. 
Then, global knowledge of communication periods and communication time instants is not 
needed as in sampled-data approaches. The use of discretized and decoupled models and the 
implementation of periodic event-triggered techniques provides a formal framework that limits 
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Fig. 5. States of four second-order agents converging to a single trajectory 
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Fig. 6. Broadcasting instants for each one of the four agents 
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Fig. 7. Broadcasting instants for each one of the four agents during the first 3 seconds of simulation 


actuation and sensing update rates and reduces communication. This method also provides the 
necessary freedom to each agent in order to determine its own broadcasting instants. 

Appendix 

Proof of Lemma [7] Sufficiency: We can express V using the following: 

N 

i =i jeATi 

Since the graph is undirected (l88l) can be written in the following form: 

V = EZi EisM I (£ T «i - &Q(i ~ (JQ(> + (JQ(i) 

= Eli Ejev, 5 (£< - (if Q((i-(i)- 

Since Q > 0 and the graph is connected it is clear that if V = 0 then = £ :/ for i,j = 1...N. 
Necessity: Consider the following expression: 

= (90) 


( 88 ) 


(89) 
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where £ represents an n —dimensional consensus state and is given by £ = by 0 where 
? = [si 9> ... SnY • Then we have: 

££ = (C 0 Q) In 0 9 = £1at 0 Qs = 0 u n- (91) 

We can conclude that if consensus is achieved then V = £ T ££ = 0. • 

Proof of Theorem [7] Since the communication graph is undirected and connected C is sym¬ 
metric and there exists a similarity transformation S such that C D — S~ 1 CS is diagonal with 
one eigenvalue equal to zero. Define T — S0I n then C D = T 1 CT = C D 0P is block diagonal 
with n eigenvalues equal to zero. 

Let us now consider the following: 

T~ l CT = T~ l [CA C + A T c C^j T 

= T~ x LTT~ X A<T + T- x A t c TT~ x LT (92) 

= C D (I N 0 A + cC D ® BF ) + (I N 0 A + cC D 0 BF) r C D . 

The term 1^ 0 A + cCd 0 BF is of the form 

A 0 
0 U 

In our case since Cd is diagonal then U is block diagonal. Furthermore, each block is given 
by Ui — A + cXiBF, i = 2...N. Then we have that (l92l) is given by the block diagonal matrix: 
diag {0 n , jC 2 , C 3 ,..., jCn}, where C, = A* [PA + A T P — 2c\ i PBB 1 P) for i — 2, ...N. Since 
c > I/A 2 and P is the solution of © we can conclude that 

\ [PA + A t P - 2c\iPBB T P) < 0, 

for % = 2...N. We can see that C has n zero eigenvalues and the rest of its eigenvalues are 
negative. 

Consider the following: 

CA c p = C0 P (I N 0 A + c£0 BF) p, (93) 
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where p is an eigenvector of £ given by p = 1 ^- < 8 > s associated with a zero eigenvalue of £ 
where ? = [51 • • • ?n] T - Then we have: 

CA c p = (£ ® PA + c£ 2 <g) PBF) p 

= C\ N ® PAe + c£ 2 1jv ® PPP^ (94) 

= OnAT- 

Similarly, AjCp = 0 nN . Then it is clear that p is an eigenvector of £ associated with a zero 
eigenvalue. • 


References 

[1] W. Ren, R. W. Beard, and E. M. Atkins, “Information consensus in multivehicle cooperative control,” IEEE Control Systems 
Magazine, vol. 27, no. 2, pp. 71-82, 2007. 

[2] Y. Cao and W. Ren, “Multi-vehicle coordination for double integrator dynamics under fixed undirected/directed interaction 
in a sampled-data setting,” International Journal of Robust and Nonlinear Control, vol. 20, pp. 987-1000, 2010. 

[3] T. Hayakawa, T. Matsuzawa, and S. Hara, “Formation control of multi-agent systems with sampled information,” in 45th 
IEEE Conference on Decision and Control, 2006, pp. 4333—4338. 

[4] H. Liu, G. Xie, and L. Wang, “Necessary and sufficient conditions for solving consensus of double integrator dynamics 
via sampled control,” International Journal of Robust and Nonlinear Control, vol. 20, no. 15, pp. 1706-1722, 2010. 

[5] J. Qin and H. Gao, “A sufficient condition for convergence of sampled-data consensus for double integrator dynamics 
with nonuniform and time-varying communication delays,” IEEE Transactions on Automatic Control, vol. 57, no. 9, pp. 
2417-2422, 2012. 

[6] K. J. Astrom, “Event based control,” in A. Astolfi and L. Marconi, (eds.), Analysis and Design of Nonlinear Control 
Systems, Springer-Verlag, Berlin, 2008, pp. 127-147. 

[7] M. C. F. Donkers and W. P. M. H. Heemels, “Output-based event-triggered control with guaranteed linfty-gain and improved 
event-triggering,” in 49th IEEE Conference on Decision and Control, 2010, pp. 3246-3251. 

[8] C. D. Persis, R. Sailer, and F. Wirth, “Parsimonious event-triggered distributed control: a zeno free approach,” Automatica, 
vol. 49, no. 7, pp. 2116-2124, 2013. 

[9] E. Garcia and P. J. Antsaklis, “Model-based event-triggered control for systems with quantization and time-varying network 
delays,” IEEE Transactions on Automatic Control, vol. 58, no. 2, pp. 422^134, 2013. 

[10] X. Wang and M. D. Lemmon, “Event-triggering in distributed networked control systems,” IEEE Transactions on Automatic 
Control, vol. 56, pp. 586-601, 2011. 

[11] E. Garcia and P. J. Antsaklis,“Parameter estimation and adaptive stabilization in time-triggered and event-triggered model- 
based control of uncertain systems,” International Journal of Control, vol. 85, no. 9, pp. 1327-1342, 2012. 

[12] E. Garcia and P. J. Antsaklis,“Output feedback model-based control of uncertain discrete-time systems with network 
induced delays,” in 51st IEEE Conference on Decision and Control, 2012, pp. 6647-6652. 

[13] P. Tabuada, “Event-triggered real-time scheduling of stabilizing control tasks,” IEEE Transactions on Automatic Control, 
vol. 52, no. 9, pp. 1680-1685, 2007. 


33 


[14] E. Garcia and P. J. Antsaklis, “Decentralized model-based event-triggered control of networked systems,” in American 
Control Conference, 2012, pp. 6485-6490. 

[15] C. Stocker, D. Vey, and J. Lunze, “Decentralized event-based control: Stability analysis and experimental evaluation,” 
Nonlinear Analysis: Hybrid Systems, vol. 10, pp. 141-155, 2013. 

[16] M. Guinaldo, D. V. Dimarogonas, K. H. Johansson, J. Sanchez, and S. Dormido, “Distributed event-based control for 
interconnected linear systems,” in 2011 50th IEEE Conference on Decision and Control and European Control Conference 
(CDC-ECC). IEEE, 2011, pp. 2553-2558. 

[17] D. V. Dimarogonas, E. Frazzoli, and K. H. Johansson, “Distributed event-triggered control for multi-agent systems,” IEEE 
Transactions on Automatic Control, vol. 57, no. 5, pp. 1291-1297, 2012. 

[18] E. Garcia, Y. Cao, H. Yu, P. J. Antsaklis, and D. W. Casbeer, “Decentralized event-triggered cooperative control with 
limited communication,” International Journal of Control, vol. 86, no. 9, pp. 1479-1488, 2013. 

[19] G. S. Seyboth, D. V. Dimarogonas, and K. H. Johansson, “Event-based broadcasting for multi-agent average consensus,” 
Automatica, vol. 49, no. 1, pp. 245-252, 2013. 

[20] H. Yu and P. J. Antsaklis, “Quantized output synchronization of networked passive systems with event-driven communi¬ 
cation,” in American Control Conference, 2012, pp. 5706-5711. 

[21] X. Yin and D. Yue, “Event-triggered tracking control for heterogeneous multi-agent systems with Markov comunication 
delays,” Journal of the Franklin Institute, vol. 350, no. 5, pp. 1312-1334, 2013. 

[22] X. Chen and F. Hao, “Event-triggered average consensus control for discrete-time multi-agent systems,” 1ET Control Theory 
and Applications, vol. 6, no. 16, pp. 2493-2498, 2012. 

[23] M. Guo and D. V. Dimarogonas, “Nonlinear consensus via continuous, sampled, and aperiodic updates,” International 
Journal of Control, vol. 86, no. 4, pp. 567-578, 2013. 

[24] A. Anta and P. Tabuada, “To sample or not to sample: Self-triggered control for nonlinear systems,” IEEE Transactions 
on Automatic Control, vol. 55, no. 9, pp. 2030-2042, 2010. 

[25] M. Mazo, A. Anta, and P. Tabuada, “An iss self-triggered implementation of linear controllers,” Automatica, vol. 46, no. 8, 
pp. 1310-1314, 2010. 

[26] C. Nowzari and J. Cortes, “Self-triggered coordination of robotic networks for optimal deployment,” Automatica, vol. 48, 
no. 6, pp. 1077-1087, 2012. 

[27] X. Wang and M. Lemmon, “State based self-triggered feedback control systems with 12 stability,” IEEE Transactions on 
Automatic Control, vol. 54, no. 3, 2009. 

[28] W. Heemels, M. Donkers, and A. R. Teel, “Periodic event-triggered control for linear systems,” Automatic Control, IEEE 
Transactions on, vol. 58, no. 4, pp. 847-861, 2013. 

[29] W. Heemels and M. Donkers, “Model-based periodic event-triggered control for linear systems,” Automatica, vol. 49, no. 3, 
pp. 698-711, 2013. 

[30] Z. Li, Z. Duan, G. Chen, and L. Huang, “Consensus of multiagent systems and synchronization of complex networks: a 
unified viewpoint,” IEEE Transactions on Circuits and Systems-I: Regular papers, vol. 57, no. 1, pp. 213-224, 2010. 

[31] Z. Li, Z. Duan, and G. Chen, “Dynamic consensus of linear multi-agent systems,” Control Theory & Applications, IET, 
vol. 5, no. 1, pp. 19-28, 2011. 

[32] C. Q. Ma and J. F. Zhang, “Necessary and sufficient conditions for consensusability of linear multi-agent systems,” IEEE 
Transactions on Automatic Control, vol. 55, no. 5, pp. 1263-1268, 2010. 


34 


[33] W. Ren, “Synchronization of coupled harmonic oscillators with local interaction,” Automatica, vol. 44, no. 12, pp. 3195— 
3200, 2008. 

[34] L. Scardovi and R. Sepulchre, “Synchronization in networks of identical linear systems,” Automatica , vol. 45, no. 11, pp. 
2557-2562, 2009. 

[35] Y. Su and J. Huang, “Stability of a class of linear switching systems with applications to two consensus problems,” IEEE 
Transactions on Automatic Control, vol. 57, no. 6, pp. 1420-1430, 2012. 

[36] S. E. Tuna, “Synchronizing linear systems via partial-state coupling,” Automatica, vol. 44. no. 8, pp. 2179-2184, 2008. 

[37] S. E. Tuna, “Conditions for synchronizability in arrays of coupled linear systems,” IEEE Transactions on Automatic Control, 
vol. 54, no. 10, pp. 2416-2420, 2009. 

[38] G. Wen, Z. Duan, W. Ren, and G. Chen, “Distributed consensus of multi-agent systems with general linear node dynamics 
and intermittent communications,” in International Journal of Robust and Nonlinear Control, 2013. 

[39] T. Liu, D. J. Hill, and B. Liu, “Synchronization of dynamical networks with distributed event-based communication,” in 
51st IEEE Conference on Decision and Control, 2012, pp. 7199-7204. 

[40] W. Zhu, Z.-R Jiang, and G. Feng, “Event-based consensus of multi-agent systems with general linear models,” Automatica, 
vol. 50, no. 2, pp. 552-558, 2014. 

[41] E. Garcia, Y. Cao, and D. W. Casbeer, “Cooperative control with general linear dynamics and limited communication: 
centralized and decentralized event-triggered control strategies,” in American Control Conference, 2014, pp. 159-164. 

[42] E. Garcia and Y. Cao and D. W. Casbeer, “Decentralized event-triggered consensus with general linear dynamics,” 
Automatica, vol. 50, no. 10, pp. 2633-2640, 2014. 

[43] O. Demir and J. Lunze, “Event-based synchronisation of multi-agent systems,” in IFAC Conference on Analysis and Design 
of Hybrid Systems, 2012, pp. 1-6. 

[44] E. Garcia, Y. Cao, and D. W. Casbeer, “Event-triggered cooperative control with general linear dynamics and communication 
delays,” in IEEE Conference on Decision and Control, 2014, pp. 2914-2919. 

[45] R. Aragues, G. Shi. D. V. Dimarogonas, C. Sagues, and K. H. Johansson, “Distributed algebraic connectivity estimation 
for adaptive event-triggered consensus,” in American Control Conference, 2012, pp. 32-37. 

[46] M. Franceschelli, A. Gasparri, A. Giua, and C. Seatzu, “Decentralized Laplacian eigenvalues estimation for networked 
multi-agent systems,” in 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference, 2009, pp. 
2717-2722. 


35 


